Coral-infecting parasites in cold marine ecosystems.

Coral reefs are a biodiversity hotspot,1,2 and the association between coral and intracellular dinoflagellates is a model for endosymbiosis.3,4 Recently, corals and related anthozoans have also been found to harbor another kind of endosymbiont, apicomplexans called corallicolids.5 Apicomplexans are a diverse lineage of obligate intracellular parasites6 that include human pathogens such as the malaria parasite, Plasmodium.7 Global environmental sequencing shows corallicolids are tightly associated with tropical and subtropical reef environments,5,8,9 where they infect diverse corals across a range of depths in many reef systems, and correlate with host mortality during bleaching events.10 All of this points to corallicolids being ecologically significant to coral reefs, but it is also possible they are even more widely distributed because most environmental sampling is biased against parasites that maintain a tight association with their hosts throughout their life cycle. We tested the global distribution of corallicolids using a more direct approach, by specifically targeting potential anthozoan host animals from cold/temperate marine waters outside the coral reef context. We found that corallicolids are in fact common in such hosts, in some cases at high frequency, and that they infect the same tissue as parasites from topical coral reefs. Parasite phylogeny suggests corallicolids move between hosts and habitats relatively frequently, but that biogeography is more conserved. Overall, these results greatly expand the range of corallicolids beyond coral reefs, suggesting they are globally distributed parasites of marine anthozoans, which also illustrates significant blind spots that result from strategies commonly used to sample microbial biodiversity.

Phylogenomics reveals Adeleorina are an ancient and distinct subgroup of Apicomplexa.

Apicomplexans are a diverse phylum of unicellular eukaryotes that share obligate relationships with terrestrial and aquatic animal hosts. Many well-studied apicomplexans are responsible for several deadly zoonotic and human diseases, most notably malaria caused by Plasmodium. Interest in the evolutionary origin of apicomplexans has also spurred recent work on other more deeply-branching lineages, especially gregarines and sister groups like squirmids and chrompodellids. But a full picture of apicomplexan evolution is still lacking several lineages, and one major, diverse lineage that is notably absent is the adeleorinids. Adeleorina apicomplexans comprises hundreds of described species that infect invertebrate and vertebrate hosts across the globe. Although historically considered coccidians, phylogenetic trees based on limited data have shown conflicting branch positions for this subgroup, leaving this question unresolved. Phylogenomic trees and large-scale analyses comparing cellular functions and metabolism between major subgroups of apicomplexans have not incorporated Adeleorina because only a handful of molecular markers and a couple organellar genomes are available, ultimately excluding this group from contributing to our understanding of apicomplexan evolution and biology. To address this gap, we have generated complete genomes from mitochondria and plastids, as well as multiple deep-coverage single-cell transcriptomes of nuclear genes from two Adeleorina species, Klossia helicina and Legerella nova, and inferred a 206-protein phylogenomic tree of Apicomplexa. We observed distinct structures reported in species descriptions as remnant host structures surrounding adeleorinid oocysts. Klossia helicina and L. nova branched, as expected, with monoxenous adeleorinids within the Adeleorina and their mitochondrial and plastid genomes exhibited similarity to published organellar adeleorinid genomes. We show with a phylogeneomic tree and subsequent phylogenomic analyses that Adeleorina are not closely related to any of the currently sampled apicomplexan subgroups, and instead fall as a sister to a large clade encompassing Coccidia, Protococcidia, Hematozoa, and Nephromycida, collectively. This resolves Adeleorina as a key independently-branching group, separate from coccidians, on the tree of Apicomplexa, which now has all known major lineages sampled.

Complete genome sequences of 12 bacterial strains from the honey bee gut, resolved with long-read nanopore sequencing.

We report the de novo sequencing of six bacterial strains isolated from the Western honey bee, as well as the resequencing of six strains that have existing draft genomes, to obtain complete, chromosomal-level assemblies. These strains include the bee gut symbiont genera Bartonella, Bifidobacterium, Snodgrassella, Gilliamella, Lactobacillus, and the opportunistic pathogen Serratia marcescens KZ11.

Reconstruction of Plastid Proteomes of Apicomplexans and Close Relatives Reveals the Major Evolutionary Outcomes of Cryptic Plastids.

Apicomplexans and related lineages comprise many obligate symbionts of animals; some of which cause notorious diseases such as malaria. They evolved from photosynthetic ancestors and transitioned into a symbiotic lifestyle several times, giving rise to species with diverse non-photosynthetic plastids. Here, we sought to reconstruct the evolution of the cryptic plastids in the apicomplexans, chrompodellids, and squirmids (ACS clade) by generating five new single-cell transcriptomes from understudied gregarine lineages, constructing a robust phylogenomic tree incorporating all ACS clade sequencing datasets available, and using these to examine in detail, the evolutionary distribution of all 162 proteins recently shown to be in the apicoplast by spatial proteomics in Toxoplasma. This expanded homology-based reconstruction of plastid proteins found in the ACS clade confirms earlier work showing convergence in the overall metabolic pathways retained once photosynthesis is lost, but also reveals differences in the degrees of plastid reduction in specific lineages. We show that the loss of the plastid genome is common and unexpectedly find many lineage- and species-specific plastid proteins, suggesting the presence of evolutionary innovations and neofunctionalizations that may confer new functional and metabolic capabilities that are yet to be discovered in these enigmatic organelles.

Host-microbiome metabolism of a plant toxin in bees.

While foraging for nectar and pollen, bees are exposed to a myriad of xenobiotics, including plant metabolites, which may exert a wide range of effects on their health. Although the bee genome encodes enzymes that help in the metabolism of xenobiotics, it has lower detoxification gene diversity than the genomes of other insects. Therefore, bees may rely on other components that shape their physiology, such as the microbiota, to degrade potentially toxic molecules. In this study, we show that amygdalin, a cyanogenic glycoside found in honey bee-pollinated almond trees, can be metabolized by both bees and members of the gut microbiota. In microbiota-deprived bees, amygdalin is degraded into prunasin, leading to prunasin accumulation in the midgut and hindgut. In microbiota-colonized bees, on the other hand, amygdalin is degraded even further, and prunasin does not accumulate in the gut, suggesting that the microbiota contribute to the full degradation of amygdalin into hydrogen cyanide. In vitro experiments demonstrated that amygdalin degradation by bee gut bacteria is strain-specific and not characteristic of a particular genus or species. We found strains of Bifidobacterium, Bombilactobacillus, and Gilliamella that can degrade amygdalin. The degradation mechanism appears to vary since only some strains produce prunasin as an intermediate. Finally, we investigated the basis of degradation in Bifidobacterium wkB204, a strain that fully degrades amygdalin. We found overexpression and secretion of several carbohydrate-degrading enzymes, including one in glycoside hydrolase family 3 (GH3). We expressed this GH3 in Escherichia coli and detected prunasin as a byproduct when cell lysates were cultured with amygdalin, supporting its contribution to amygdalin degradation. These findings demonstrate that both host and microbiota can act together to metabolize dietary plant metabolites.

Most plants produce chemicals that are toxic to at least some animals. Whether or not the toxins are harmful to a particular animal depends on how much they consume and the specific biochemistry that occurs during digestion. The enzymes produced in the gut both by the animal and by the microbes that reside there often help break down toxic substances into less harmful molecules. However, some products of this breakdown can be toxic themselves. While these products can harm the animal, they may also be detrimental to parasites living in the gut, resulting in an overall positive effect. Almonds and their pollen are consumed by humans and bees without apparent harmful effects. However, almonds contain amygdalin, a molecule that can produce the highly toxic compound hydrogen cyanide upon digestion. Although amygdalin can be toxic to bees in high doses, the amount usually found in almond nectar is not harmful, and indeed, it may protect bees from parasites. Motta et al. wanted to know how amygdalin is digested in the gut of bees, and whether gut microbes have a role in this digestion. To answer these questions, Motta et al. compared the effects of consuming amygdalin on normal bees and bees lacking gut microbes. Bees without gut microbes broke down amygdalin into a harmless substance called prunasin. However, only bees with gut microbes could further break down prunasin into hydrogen cyanide. Interestingly, the full metabolism of amygdalin had no detectable effect on whether the bees survived for longer times or on which microbes were found in the gut. Motta et al. also found some gut bacteria in bees that can break down amygdalin and release hydrogen cyanide, and identified the enzyme responsible for the process. When the gene encoding this enzyme was inserted into a different species of bacteria, the second species gained the ability to break down amygdalin. The findings of Motta et al. explain a role of gut microbes in processing amygdalin in bees. In the future, this may be the key to understanding how humans and other creatures process plant toxins. Future work on the relationship between animals and microbes living in their guts could help scientists understand how to manipulate the digestion and processing of toxins, nutrients, or drugs to benefit human health.

Gene Transfer Agents in Bacterial Endosymbionts of Microbial Eukaryotes.

Gene transfer agents (GTAs) are virus-like structures that package and transfer prokaryotic DNA from donor to recipient prokaryotic cells. Here, we describe widespread GTA gene clusters in the highly reduced genomes of bacterial endosymbionts from microbial eukaryotes (protists). Homologs of the GTA capsid and portal complexes were initially found to be present in several highly reduced alphaproteobacterial endosymbionts of diplonemid protists (Rickettsiales and Rhodospirillales). Evidence of GTA expression was found in polyA-enriched metatranscriptomes of the diplonemid hosts and their endosymbionts, but due to biases in the polyA-enrichment methods, levels of GTA expression could not be determined. Examining the genomes of closely related bacteria revealed that the pattern of retained GTA head/capsid complexes with missing tail components was common across Rickettsiales and Holosporaceae (Rhodospirillales), all obligate symbionts with a wide variety of eukaryotic hosts. A dN/dS analysis of Rickettsiales and Holosporaceae symbionts revealed that purifying selection is likely the main driver of GTA evolution in symbionts, suggesting they remain functional, but the ecological function of GTAs in bacterial symbionts is unknown. In particular, it is unclear how increasing horizontal gene transfer in small, largely clonal endosymbiont populations can explain GTA retention, and, therefore, the structures may have been repurposed in endosymbionts for host interactions. Either way, their widespread retention and conservation in endosymbionts of diverse eukaryotes suggests an important role in symbiosis.

Taxonomy of the Apicomplexan Symbionts of Coral, including Corallicolida ord. nov., Reassignment of the Genus Gemmocystis, and Description of New Species Corallicola aquarius gen. nov. sp. nov. and Anthozoaphila gnarlus gen. nov. sp. nov.

Corals (Metazoa; Cnidaria; Anthozoa) have recently been shown to play host to a widespread and diverse group of intracellular symbionts of the phylum Apicomplexa. These symbionts, colloquially called "corallicolids", are mostly known through molecular analyses, and no formal taxonomy has been proposed. Another apicomplexan, Gemmocystis cylindrus (described from the coral Dendrogyra cylindrus), may be related to corallicolids, but lacks molecular data. Here, we isolate and describe motile trophozoite (feeding) corallicolids cells using microscopic (light, SEM, and TEM) and molecular phylogenetic analysis to provide the basis for a formal description. Phylogenetic analyses using nuclear and plastid rRNA operons, and three mitochondrial protein sequences derived from single cell transcriptomes, all confirm that these organisms fall into a discrete deep-branching clade within the Apicomplexa not closely related to any known species or major subgroup. As a result, we assign this clade to a new order, Corallicolida ord. nov., and family, Corallicolidae fam. nov. We describe a type species, Corallicola aquarius gen. nov. sp. nov. from its Rhodactis sp. host, and also describe a second species, Anthozoaphila gnarlus gen. nov. sp. nov., from the coral host Madracis mirabilis. Finally, we propose reassigning the incertae sedis taxon G. cylindrus from the order Agamococcidiorida to the Corallicolida, based on similarities in morphology and host localization to that of the corallicolids.

Phylogenomics Identifies a New Major Subgroup of Apicomplexans, Marosporida class nov., with Extreme Apicoplast Genome Reduction.

The phylum Apicomplexa consists largely of obligate animal parasites that include the causative agents of human diseases such as malaria. Apicomplexans have also emerged as models to study the evolution of nonphotosynthetic plastids, as they contain a relict chloroplast known as the apicoplast. The apicoplast offers important clues into how apicomplexan parasites evolved from free-living ancestors and can provide insights into reductive organelle evolution. Here, we sequenced the transcriptomes and apicoplast genomes of three deep-branching apicomplexans, Margolisiella islandica, Aggregata octopiana, and Merocystis kathae. Phylogenomic analyses show that these taxa, together with Rhytidocystis, form a new lineage of apicomplexans that is sister to the Coccidia and Hematozoa (the lineages including most medically significant taxa). Members of this clade retain plastid genomes and the canonical apicomplexan plastid metabolism. However, the apicoplast genomes of Margolisiella and Rhytidocystis are the most reduced of any apicoplast, are extremely GC-poor, and have even lost genes for the canonical plastidial RNA polymerase. This new lineage of apicomplexans, for which we propose the class Marosporida class nov., occupies a key intermediate position in the apicomplexan phylogeny, and adds a new complexity to the models of stepwise reductive evolution of genome structure and organelle function in these parasites.

Microbiome Evolution: Having the Guts to Be Different.

Metagenomic sequencing of the gut microbial communities of two closely related bee species, the Western honey bee (Apis mellifera) and the Eastern honey bee (Apis cerana), show that organisms with similar characteristics can harbor unexpected differences in their microbiomes.

Highly Reduced Genomes of Protist Endosymbionts Show Evolutionary Convergence.

Genome evolution in bacterial endosymbionts is notoriously extreme: the combined effects of strong genetic drift and unique selective pressures result in highly reduced genomes with distinctive adaptations to hosts [1-4]. These processes are mostly known from animal endosymbionts, where nutritional endosymbioses represent the best-studied systems. However, eukaryotic microbes, or protists, also harbor diverse bacterial endosymbionts, but their genome reduction and functional relationships with their hosts are largely unexplored [5-7]. We sequenced the genomes of four bacterial endosymbionts from three species of diplonemids, poorly studied but abundant and diverse heterotrophic protists [8-12]. The endosymbionts come from two bacterial families, Rickettsiaceae and Holosporaceae, that have invaded two families of diplonemids, and their genomes have converged on an extremely small size (605-632 kilobase pairs [kbp]), similar gene content (e.g., metabolite transporters and secretion systems), and reduced metabolic potential (e.g., loss of energy metabolism). These characteristics are generally found in both families, but the diplonemid endosymbionts have evolved greater extremes in parallel. They possess modified type VI secretion systems that could function in manipulating host metabolism or other intracellular interactions. Finally, modified cellular machinery like the ATP synthase without oxidative phosphorylation, and the reduced flagellar apparatus present in some diplonemid endosymbionts and nutritional animal endosymbionts, indicates that intracellular mechanisms have converged in bacterial endosymbionts with various functions and from different eukaryotic hosts across the tree of life.

Division of labor in honey bee gut microbiota for plant polysaccharide digestion.

Bees acquire carbohydrates from nectar and lipids; and amino acids from pollen, which also contains polysaccharides including cellulose, hemicellulose, and pectin. These potential energy sources could be degraded and fermented through microbial enzymatic activity, resulting in short chain fatty acids available to hosts. However, the contributions of individual microbiota members to polysaccharide digestion have remained unclear. Through analysis of bacterial isolate genomes and a metagenome of the honey bee gut microbiota, we identify that Bifidobacterium and Gilliamella are the principal degraders of hemicellulose and pectin. Both Bifidobacterium and Gilliamella show extensive strain-level diversity in gene repertoires linked to polysaccharide digestion. Strains from honey bees possess more such genes than strains from bumble bees. In Bifidobacterium, genes encoding carbohydrate-active enzymes are colocated within loci devoted to polysaccharide utilization, as in Bacteroides from the human gut. Carbohydrate-active enzyme-encoding gene expressions are up-regulated in response to particular hemicelluloses both in vitro and in vivo. Metabolomic analyses document that bees experimentally colonized by different strains generate distinctive gut metabolomic profiles, with enrichment for specific monosaccharides, corresponding to predictions from genomic data. The other 3 core gut species clusters (Snodgrassella and 2 Lactobacillus clusters) possess few or no genes for polysaccharide digestion. Together, these findings indicate that strain composition within individual hosts determines the metabolic capabilities and potentially affects host nutrition. Furthermore, the niche specialization revealed by our study may promote overall community stability in the gut microbiomes of bees.

A widespread coral-infecting apicomplexan with chlorophyll biosynthesis genes.

Apicomplexa is a group of obligate intracellular parasites that includes the causative agents of human diseases such as malaria and toxoplasmosis. Apicomplexans evolved from free-living phototrophic ancestors, but how this transition to parasitism occurred remains unknown. One potential clue lies in coral reefs, of which environmental DNA surveys have uncovered several lineages of uncharacterized basally branching apicomplexans1,2. Reef-building corals have a well-studied symbiotic relationship with photosynthetic Symbiodiniaceae dinoflagellates (for example, Symbiodinium3), but the identification of other key microbial symbionts of corals has proven to be challenging4,5. Here we use community surveys, genomics and microscopy analyses to identify an apicomplexan lineage-which we informally name 'corallicolids'-that was found at a high prevalence (over 80% of samples, 70% of genera) across all major groups of corals. Corallicolids were the second most abundant coral-associated microeukaryotes after the Symbiodiniaceae, and are therefore core members of the coral microbiome. In situ fluorescence and electron microscopy confirmed that corallicolids live intracellularly within the tissues of the coral gastric cavity, and that they possess apicomplexan ultrastructural features. We sequenced the genome of the corallicolid plastid, which lacked all genes for photosystem proteins; this indicates that corallicolids probably contain a non-photosynthetic plastid (an apicoplast6). However, the corallicolid plastid differs from all other known apicoplasts because it retains the four ancestral genes that are involved in chlorophyll biosynthesis. Corallicolids thus share characteristics with both their parasitic and their free-living relatives, which suggests that they are evolutionary intermediates and implies the existence of a unique biochemistry during the transition from phototrophy to parasitism.

Is Host Filtering the Main Driver of Phylosymbiosis across the Tree of Life?

Host-associated microbiota composition can be conserved over evolutionary time scales. Indeed, closely related species often host similar microbiota; i.e., the composition of their microbiota harbors a phylogenetic signal, a pattern sometimes referred to as "phylosymbiosis." Elucidating the origins of this pattern is important to better understand microbiota ecology and evolution. However, this is hampered by our lack of theoretical expectations and a comprehensive overview of phylosymbiosis prevalence in nature. Here, we use simulations to provide a simple expectation for when we should expect this pattern to occur and then review the literature to document the prevalence and strength of phylosymbiosis across the host tree of life. We demonstrate that phylosymbiosis can readily emerge from a simple ecological filtering process, whereby a given host trait (e.g., gut pH) that varies with host phylogeny (i.e., harbors a phylogenetic signal) filters preadapted microbes. We found marked differences between methods used to detect phylosymbiosis, so we proposed a series of practical recommendations based on using multiple best-performing approaches. Importantly, we found that, while the prevalence of phylosymbiosis is mixed in nature, it appears to be stronger for microbiotas living in internal host compartments (e.g., the gut) than those living in external compartments (e.g., the rhizosphere). We show that phylosymbiosis can theoretically emerge without any intimate, long-term coevolutionary mechanisms and that most phylosymbiosis patterns observed in nature are compatible with a simple ecological process. Deviations from baseline ecological expectations might be used to further explore more complex hypotheses, such as codiversification. IMPORTANCE Phylosymbiosis is a pattern defined as the tendency of closely related species to host microbiota whose compositions resemble each other more than host species drawn at random from the same tree. Understanding the mechanisms behind phylosymbiosis is important because it can shed light on rules governing the assembly of host-associated microbiotas and, potentially, their coevolutionary dynamics with hosts. For example, is phylosymbiosis a result of coevolution, or can it be generated by simple ecological filtering processes? Beyond qualitative theoretical models, quantitative theoretical expectations can provide new insights. For example, deviations from a simple baseline of ecological filtering may be used to test more-complex hypotheses (e.g., coevolution). Here, we use simulations to provide evidence that simple host-related ecological filtering can readily generate phylosymbiosis, and we contrast these predictions with real-world data. We find that while phylosymbiosis is widespread in nature, phylosymbiosis patterns are compatible with a simple ecological model in the majority of taxa. Internal compartments of hosts, such as the animal gut, often display stronger phylosymbiosis than expected from a purely ecological filtering process, suggesting that other mechanisms are also involved.

Competitive organelle-specific adaptors recruit Vps13 to membrane contact sites.

The regulated expansion of membrane contact sites, which mediate the nonvesicular exchange of lipids between organelles, requires the recruitment of additional contact site proteins. Yeast Vps13 dynamically localizes to membrane contacts that connect the ER, mitochondria, endosomes, and vacuoles and is recruited to the prospore membrane in meiosis, but its targeting mechanism is unclear. In this study, we identify the sorting nexin Ypt35 as a novel adaptor that recruits Vps13 to endosomal and vacuolar membranes. We characterize an interaction motif in the Ypt35 N terminus and identify related motifs in the prospore membrane adaptor Spo71 and the mitochondrial membrane protein Mcp1. We find that Mcp1 is a mitochondrial adaptor for Vps13, and the Vps13-Mcp1 interaction, but not Ypt35, is required when ER-mitochondria contacts are lost. All three adaptors compete for binding to a conserved six-repeat region of Vps13 implicated in human disease. Our results support a competition-based model for regulating Vps13 localization at cellular membranes.

Author Correction: Convergent evolution of a modified, acetate-driven TCA cycle in bacteria.

In this Brief Communication, the authors omitted references to several previous studies that have demonstrated that the TCA variant shown has also been found in several other bacterial species, specifically among some anaerobic Deltaproteobacteria. The Brief Communication focused on showing that this variant is more widespread than previously known, particularly in animal-associated bacteria. Using genetic approaches, Kwong et al. demonstrated this alternative cycle in additional, distantly related organisms. Thus, the previous work does not affect the results of the advance provided. However, the previous studies are relevant to this topic and should have been cited, and Kwong et al. apologize for this for the omission. Citations to several of the studies should have been included as refs 9-13 in the text as shown.

Genome Sequences of Apibacter spp., Gut Symbionts of Asian Honey Bees.

Honey bees have distinct gut microbiomes consisting almost entirely of several host-specific bacterial species. We present the genomes of three strains of Apibacter spp., bacteria of the Bacteroidetes phylum that are endemic to Asian honey bee species (Apis dorsata and Apis cerana). The Apibacter strains have similar metabolic abilities to each other and to Apibacter mensalis, a species isolated from a bumble bee. They use microaerobic respiration and fermentation to catabolize a limited set of monosaccharides and dicarboxylic acids. All strains are capable of gliding motility and encode a type IX secretion system. Two strains and A. mensalis have type VI secretion systems, and all strains encode Rhs or VgrG proteins used in intercellular interactions. The characteristics of Apibacter spp. are consistent with adaptions to life in a gut environment; however, the factors responsible for host-specificity and mutualistic interactions remain to be uncovered.

Microbiome Structure Influences Infection by the Parasite Crithidia bombi in Bumble Bees.

Recent declines in bumble bee populations are of great concern and have prompted critical evaluations of the role of pathogen introductions and host resistance in bee health. One factor that may influence host resilience when facing infection is the gut microbiota. Previous experiments with Bombus terrestris, a European bumble bee, showed that the gut microbiota can protect against Crithidia bombi, a widespread trypanosomatid parasite of bumble bees. However, the particular characteristics of the microbiome responsible for this protective effect have thus far eluded identification. Using wild and commercially sourced Bombus impatiens, an important North American pollinator, we conducted cross-wise microbiota transplants to naive hosts of both backgrounds and challenged them with a Crithidia parasite. As with B. terrestris, we find that microbiota-dependent protection against Crithidia operates in B. impatiens Lower Crithidia infection loads were experimentally associated with high microbiome diversity, large gut bacterial populations, and the presence of Apibacter, Lactobacillus Firm-5, and Gilliamella spp. in the gut community. These results indicate that even subtle differences between gut community structures can have a significant impact on a microbiome's ability to defend against parasite infections.IMPORTANCE Many wild bumble bee populations are under threat due to human activity, including through the introduction of pathogens via commercially raised bees. Recently, it was found that the bumble bee gut microbiota can help defend against a common parasite, Crithidia bombi, but the particular factors contributing to this protection are unknown. Using both wild and commercially raised bees, we conducted microbiota transplants to show that microbiome diversity, total gut bacterial load, and the presence of certain core members of the microbiota may all impact bee susceptibility to Crithidia infection. Bee origin (genetic background) was also a factor. Finally, by examining this phenomenon in a previously uninvestigated bee species, our study demonstrates that microbiome-mediated resistance to Crithidia is conserved across multiple bumble bee species. These findings highlight how intricate interactions between hosts, microbiomes, and parasites can have wide-ranging consequences for the health of ecologically important species.

Diversification of Type VI Secretion System Toxins Reveals Ancient Antagonism among Bee Gut Microbes.

Microbial communities are shaped by interactions among their constituent members. Some Gram-negative bacteria employ type VI secretion systems (T6SSs) to inject protein toxins into neighboring cells. These interactions have been theorized to affect the composition of host-associated microbiomes, but the role of T6SSs in the evolution of gut communities is not well understood. We report the discovery of two T6SSs and numerous T6SS-associated Rhs toxins within the gut bacteria of honey bees and bumble bees. We sequenced the genomes of 28 strains of Snodgrassella alvi, a characteristic bee gut microbe, and found tremendous variability in their Rhs toxin complements: altogether, these strains appear to encode hundreds of unique toxins. Some toxins are shared with Gilliamella apicola, a coresident gut symbiont, implicating horizontal gene transfer as a source of toxin diversity in the bee gut. We use data from a transposon mutagenesis screen to identify toxins with antibacterial function in the bee gut and validate the function and specificity of a subset of these toxin and immunity genes in Escherichia coli Using transcriptome sequencing, we demonstrate that S. alvi T6SSs and associated toxins are upregulated in the gut environment. We find that S. alvi Rhs loci have a conserved architecture, consistent with the C-terminal displacement model of toxin diversification, with Rhs toxins, toxin fragments, and cognate immunity genes that are expressed and confer strong fitness effects in vivo Our findings of T6SS activity and Rhs toxin diversity suggest that T6SS-mediated competition may be an important driver of coevolution within the bee gut microbiota.IMPORTANCE The structure and composition of host-associated bacterial communities are of broad interest, because these communities affect host health. Bees have a simple, conserved gut microbiota, which provides an opportunity to explore interactions between species that have coevolved within their host over millions of years. This study examined the role of type VI secretion systems (T6SSs)-protein complexes used to deliver toxic proteins into bacterial competitors-within the bee gut microbiota. We identified two T6SSs and diverse T6SS-associated toxins in bacterial strains from bees. Expression of these genes is increased in bacteria in the bee gut, and toxin and immunity genes demonstrate antibacterial and protective functions, respectively, when expressed in Escherichia coli Our results suggest that coevolution among bacterial species in the bee gut has favored toxin diversification and maintenance of T6SS machinery, and demonstrate the importance of antagonistic interactions within host-associated microbial communities.